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(1.1) 



Abstract. We consider the Wigner-Fokker-Planck equation subject to con- 
fining potentials which are smooth perturbations of the harmonic oscillator. 
For a certain class of these perturbations we prove that the equation admits 
a unique stationary solution in a weighted L 2 -space. Moreover, we show that 
the solutions of the time-dependent problem converge towards this steady state 
with an exponential rate. 



1. Introduction 

This work is devoted to the study of the Wigner-Fokker-Planck equation (WFP), 
considered in the following dimensionless form 

f d t w + £ ■ V x w + 0[V]w = Ai=w + 2 div c (£_w) + A x w, 

\ iu\ t=0 = wo(x,Q, 

where x,£ G R d , for d > 1, and t G K + . Here, w(t, x, £) is the Wigner transform 
[22j of a quantum mechanical density matrix p(t, x, y), as defined by 

(1.2) ^> : =(^/ R /(^ + I >*-I)e^. 

Recall that, for any time t G K+, a quantum mechanical (mixed) state is given by 
a positive, self-adjoint trace class operator p(t) G 3^~, where 

Ti := {p G Q3(L 2 (M d )) : tr \p\ < oo} . 

Since C T2, the space of Hilbert- Schmidt operators, i.e. 

T 2 := {/O S < B(L 2 (M. d )) : ti(p*p) < 00} , 

we can identify p(t) with its corresponding integral kernel pit,-,-) E L 2 (R 2d ), the 
so-called density matrix. Consequently, p{t) acts on a function tp £ i 2 (R d ) via 



(p(t)<p)(x)= p(t,x,y)ip{y)dy. 



Using the Wigner transformed density matrix (|1.2|) , one obtains a phase-space de- 
scription of quantum mechanics, with x being the position and £ the momentum. 
By definition, w(t, -, •) G i 2 (M 2d ) and real-valued. In contrast to classical phase 
space distributions, w(t, x, £) in general also takes negative values. 

Equation (|l.ip governs the time evolution of w(t,x,£) in the framework of so- 
called open quantum systems, which model both the Hamiltonian evolution of a 
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quantum system and its interaction with an environment (see [12] . e.g.). Here, we 
specifically describe these interactions by the Fokker-Planck type diffusion operator 
on the r.h.s of (jl.ip . Potential forces acting on w(t, ■, •) are taken into account by 
the pseudo-differential operator 

(1.3) (®[V]f)(x,Q := — ^ ff SV{x,r,)f(x,t') e^"^ d£' dr,, 
where the symbol SV is given by 

(1.4) sV(x,r,) = v(x + ^)-v(x-^j. 

The WFP equation is a kinetic model for quantum mechanical charge-transport, 
including diffusive effects, as needed, e.g., in the description of quantum Brownian 
motion [14j . quantum optics [16j . and semiconductor device simulations [15j . It 
can be considered as a quantum mechanical generalization of the usual kinetic 
Fokker-Planck equation (or Kramer's equation), to which it is known to converge 
in the classical limit h — > 0, after an appropriate rescaling of the appearing physical 
parameters [§]. The WFP has been partly derived in [TO], as a rigorous scaling 
limit for a system of particles interacting with a heat bath of phonons. Formal 
derivations (based on physical arguments) can also can be found in [T31 [TU, H?T] . 

In recent years, mathematical studies mainly focused on the Cauchy problem of 
WFP type equations -both linear and nonlinear- in various functional analytical 
settings, see G2 HI El IH1 [UJ ■ To this end it is important to note that the dynamics 
induced by (|1.1[) maps if (L 2 (R d )) into itself, since the so-called Lindblad condition 
is fulfilled (see Remark ll. 31 below). For more details on this we refer to Eld] and the 
references given therein. In the present work we shall mainly be be concerned with 
the steady state problem corresponding to . Stationary problems for quantum 
systems based on the Wigner formalism seem to be rather difficult, as only very 
few results exist (in spite of significant efforts) : The stationary, inflow-problem for 
the linear Wigner equation in d = 1 was analyzed in [5]. The only existing result 
for the WFP equation is [5Dj, where the existence of a unique steady state for a 
quadratic potential V has been proved. We remark that criteria for the existence of 
stationary density matrices for quantum dynamical semi-groups (corresponding to 
a Heisenberg-von Neumann evolution equation in Lindblad form) were given in [17] . 
The applicability of those results to the WFP equation, however, is not obvious and 
has not yet been investigated. 

More specifically, we shall from now on assume that the potential V, appearing 
in (II. 1|) , is of the form 

(1.5) V(x) = ^\x\ 2 + XVo(x), 

with A £ M some given perturbation constant. In other words we consider a per- 
turbation Vo of the harmonic oscillator. An easy calculation shows that for such a 
V the stationary equation, corresponding to (|1.1[) , can be written as 

(1.6) Lw = \G[V ]w, 
where L is the linear operator 

(1.7) Lw := -£ ■ X7 x w + x ■ V c w + A ( w + 2div 4 (^w) + A x w. 

To solve (|1.6j) we shall use a fixed point iteration for w. However, L has a non- 
trivial kernel. Indeed it is been proved in [5U], that in the case A = there exists a 
unique stationary solution /i G 6>(R 2d ), satisfying 



(1.8) 



Lfi = 0, 
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together with the normalization condition 

(1.9) ff ii(x,Z)dxd£ = l. 

Remark 1.1. Note that for any p S Ti such that w E L 1 (M 2d ) the following formal 
identity 

trp= / p(x,x)dx — // w(x, £) dxd£, 
Jm d J J«_ 2d 

can be rigorously justified by a limiting procedure in Ti, see pQ. Since tip is 

proportional to the total mass of the quantum system, we can interpret condition 

(|1.9p as a mass normalization. 

Explicitly, p can be written as 

p = e~ A ^\ 

where the function A is given by 

A(x,0 : =I(|af + 2x.£ + 3|£| 2 )+c, 

and the constant u> E K is chosen such that JJ R2 d p(x, £,) dx d£ = 1 holds. In the 
following, we denote by a > the biggest constant such that 

(1.10) HcssA(:c,£) -crl > 0, for all (x,£) E R 2d , 

in the sense of positive definite matrices, where I is the identity matrix on M. 2d . The 
constant a > can be explicitly computed in our case to be equal to a = 1 — l/y/2. 

Next, let us define the weighted Hilbert space Ti — L 2 (R 2d , p^ 1 dxd£) equipped 
with the inner product 

(f,9) H = ff — dxdC. 

J JVL 2d M 

In the following, we also denote by 

<p(k) ~ f ip(x)e- ik - x dx, 

JR d 

the Fourier transform with respect to the variable x. With these definitions at 
hand, we can now state the first main result of our work. 

Theorem 1. Let V (x) E C°°(R d ) such that 

d ,tr, 11 v 2 ) L 'i"> 

Also, let \X\ < a /Cq, where a > is defined in (|1 . lOj) . Then, it holds: 

(i) Equation (jl.6p admits a unique weak solution Woo E Ti PI i/ 1 (M 2d ), satis- 
fying the normalization condition J f R2d Woo dx d£ = 1 . Moreover, we have 
Woc EH 2 oc (R 2d ). 

(ii) For any initial data Wq E Ti of satisfying JJ R2 d Wq dxdt; = 1, we have 

II w(t) - Woo \\ H < e~ Kt || w - Woo \\ n , t > 0, 

wftere k := ^2(C |A| - a) > 0. 

Roughly speaking, Theorem 1 shows existence and exponential stability of the 
steady state Woo for smooth perturbations Vb, which are not "too concentrated" for 
x E R d . 
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Remark 1.2. A possible class of admissible potentials Vo, satisfying condition 
(jl.lip , is given by a convolution of the form 

where, as long as e > d, one may choose any U G C(R d ), such that U € L 1 (M d ). 

Theorem[T]is formulated in the Wigner transformed picture of quantum mechan- 
ics. We shall now turn our attention to the corresponding density matrix operators 
p(t). This is important since it is a-priori not clear that Woo has a physical mean- 
ingful interpretation in terms of a positive trace class operator. To this end we 
denote by poo the Hilbcrt-Schmidt operator corresponding to the kernel p 00 (x,y), 
which is obtained from w oc (x,(,) by the inverse Wigner transform, i.e. 

Poo{x,y) = J Wee «' * " " d£. 

Analogously we denote by po the Hilbert-Schmidt operator corresponding to wq £ 
H, the initial Wigner function. 

Theorem 2. Under the same assumptions on Vq and wq as before, we have: 

(i) The steady state p^ is a positive trace-class operator on L 2 (M. d ), satisfying 
trpoo = 1. 

(ii) The steady state p^ is exponentially stable, in the sense that 

II p(t) - Poo || Xa < i^c- Kt || w - Woo \\ H , t > 0, 

for some K > depending only on the spatial dimension d. 

(iii) // the initial state u>o € Ti corresponds to a density matrix po € ( and 
hence trpo = // wo da;d^ = 1), then 

lim || p{t) - pooll^ = 0. 

Note that - in the presented framework - we do not obtain exponential conver- 
gence towards the steady state in the Ti-norm but only in the sense of Hilbert- 
Schmidt operators. 

Remark 1.3. It is straightforward to extend our results to the more general Fokker- 
Planck type operator appearing on the r.h.s. of (jl.ip . 

Qui := aA^w + 2/3div a; (V 5 w) + 27div 5 (£w) + SA x w, 

as long as the Lindblad condition holds, i.e. 

-v 2 

(1.12) a>0, ad - i/3 2 + -j) > 0. 

This implies that discarding in the diffusion in x, and hence reducing the 

r.h.s. to the classical Fokker-Planck operator Q cl w :— A^w + 2div^(^w), would not 
satisfy (|1.12|) . Nevertheless, this is a frequently used model in applications [23], 
yielding reasonable results. 

The paper is organized as follows: In Section [2] we first prove several preliminary 
results on the unperturbed problem to be used later on. In Section [3] we first 
estimate the perturbation operator 0[Vo] and then give the proofs of Theorem [1] 
and Theorem [2J 
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2. Preliminary results on the unperturbed problem 

It has been shown in (2U] that the operator L, defined in (|1.7jl . can be rewritten 
in the following form 

(2.1) Lw = div (Vw + w(V A + F)) , 

where from now on all differential operators act with respect to both x and £ (if 
not indicated otherwise). In (|2.1[) , the function A is defined by (JTJ) and 

F: =(*+k)- vA 



The reason to do so is that (|2.ip falls into the class of non-symmetric Fokker- 
Planck type operators considered in [7]. To proceed further we decompose L into 
its symmetric and anti-symmetric part in TL, i.e. we write 

(2.2) L = L S + L as 1 
where 

L s w = div (Vw + wWA) , L as w = div(Fw). 
In particular, it has been shown in [20j . that the following property holds: 

(2.3) dbr(Ffj) = 0. 

Next we shall properly define the operators L and L* . To this end we first consider 
L|„„, which is closable (w.r.t. the 7i-norm) due to its dissipativity (cf. Lemma l2~2l 

below). The operator L\ is now closed and densely defined on TL. The following 



lemma then states an important assertion for this procedure. 

Lemma 2.1. Let P = p(x, £, V x , Vg), where p is a second order polynomial. Then 

P „„ is the maximum extension of P inTL. 
io D 

Proof. The proof is generalization of the one given in [5], where the analogous 
lemma was proved on L 2 (M. 2d ). Generalizations to L 2 -spaces with polynomial 
weights in £ can be found in [H [3] . 

For all / G TL with Pf € 7i, we need to construct an approximating sequence 
{/„} ll£ N € C °°(R 2<i ), such that /„ -> f in the graph norm ||/|| P := ||/||„ + ||P/|| W . 
To this end we rewrite the weight in H as 

^{x, = e«£ ^ (N 2 + 2x • { + 3|^|T , 

since all terms appearing in this expansion can then be treated analogously to those 
in [3]. The proof however is rather awkward and lengthy and we shall leave the 
details to the reader. □ 



Since L has this polynomial form, we define L :— L „„ , obtaining the maximally 
extended L on D(L) = {/ E TL : Lf £ 7i}. For this operator L we have the following 
result. 

Lemma 2.2. L is dissipative, i.e. it holds (Lw,w) n < 0, for all w £ F){L). 
Proof. Taking into account that VA = — /i _1 V/u we have, on the one hand 

(L s w,w) H = // — div (Vw + wWA) dx d£ 
"div^V (^Y)dxd* 



2d /i \ \ /X 

2 







/ * 









dx d£ < 0. 
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On the other hand, from (|2.3p it follows 



and thus 



w divF = F ■ V/j 



div(Fw) = —fiF ■ [ 4V/i - —) — ■ V (— 
\M 2 M / \V 



An easy calculation then shows that 

(L as w,w) n = ff -div(Fw)dxd£= ff -F- V (- ] fidxdt, 
J jR 2d M J JR 2d M \M/ 

by (|2.3[) . In summary we have shown that {Lw,w) H < holds. □ 

Next we consider L*| cgo , defined by = (f,L*g) n , for /, 5 e C °°(M 2d ). 

One explicitly finds that L* — L s — i as , on C5°(lR 2d ), and thus we also have that 
L* is dissipative. Since Lemma 12.11 likewise applies to L* := L*\ we infer 
that L* is dissipative on all of its domain. 

Corollary 2.3. L generates a Co-semigroup of contractions on Ti. 

Proof. Since L and L* are both dissipative operators, it follows from the Lumer- 
Phillips theorem that L is the generator of a Co-semigroup of contractions on Ti, 
see Corollary 1.4.4 in [15]. □ 

Let, for any t > 0, e Lt be the Co-semigroup on Ti generated by L. Obviously 
we have e Lt /i = //, by ()1.8|) . In [20] the Green's function of e Lt was computed 
explicitly It shows that e Lt conserves mass, i.e. for w(t) = e Lt wo it holds 

// w(t,x,£)dxd£ = // w (x,£)dxd£, Vt > 0. 

We remark though that the functional analytical setting in [5D] is not the same as in 
the present work and thus we can not directly conclude from [5D] , that L generates 
a Co-semigroup on H. 
To proceed we define 

n 1 - := {/ e H : / i. n}, 

which is a closed subset of Ti. Note that / _L (i simply means that 
{f,H) H = ff f(x 1 0dxd£ = Q. 

J JR 2d 

It is then easily seen that e Lt maps Ti' 1 into itself: For fo e Ti.- 1 the conservation 
of mass implies 

(e Lt f Q ,vL) n = ff f(t,x,t)dxd£= ff f (x,0dxd£ = 0, Vt>0. 

J jR 2d J jR 2d 

Now, Corollary 12 . 31 allows us to prove the following key result. 
Proposition 2.4. It holds 

\\L 1 \\b(h j -) — ~ ■ 



where a > is defined in (|1 .10[) . 
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Proof. Condition (jl.lOp implies, c.f. [7, Theorem 3.10], that L has a spectral gap 
of size a > 0, and hence 



|| e^iuo - Woo)\\ H < e CTt || w - mIIw 

Here, u>o satisfies JJ R2d wo dxd£ = JJ R2d ^ dxd£ = 1. The Hille-Yosida theorem 
applied to e Lt \ H± consequently yields 



\\(L-z) 1 \\ B (H±) ^ 7X77' Vz>-ct, 



z + a 

which proves the assertion for z = 0. 

3. Proof of the main results 



□ 



Before we turn to the proof of our main theorems we need to derive an estimate 
on the perturbation induced by V$. The following lemma helps us to do so. 

Lemma 3.1. For any £ € R d we can estimate 

d 

| C | 2 „ < jn-x J2 VneN. 

m=l 

Proof. Denote by 3 := {&}£, =1 and || 3 \\% := \Y? m =i IU| 2p , for p G N. Then, 
using Holder's inequality we obtain 

iei 2 - = d»ii 3 117, < d-w 3 n» iiin?,,^,) = d- 1 )f e , 

m— 1 

which is the assertion of the lemma. □ 

With the above given lemma in hand we shall now derive an appropriate bound 
on 0[Vb] in H. 

Proposition 3.2. Let Vq(x) € C°°(M d ) be such that condition p. lip holds. Then 
the operator Q[Vq] maps Tt into Tt and it holds 

II Q[V ] \\ B(H) < C . 
Proof. First we note that /i -1 < e'^' + i\ x \ and hence 

mV }w\\ n <\\(e[Vo}w)e^ 2+ ^ 2 \\ L2{m2d) . 
In the following we denote by 

{T^ v w){x, rf) = w{x, n) := / w(x, f) e~ ie ''d£, 

the partial Fourier transform with respect to the variable (el 1 ' only and recall 
from (|1.3|) . that the operator 0[Vb] acts via 

Q[V Q ]w = T-^(SV (F^ v w))- 

Setting u(x, £) := e^l^l w(x, £) and using the series expansion of the exponential, 
we obtain 



\\(e[v }w)e^ 2+ i^ 2 \\ L2{R2d) = \\(e[v }u)e^ 



(e[v ]u) 



L 2 (R 2d 



n! 

riGNo 



L 2 (R 2d ) 
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Now, we can apply Lemma [3~T1 to estimate 

l€l a " 



(3.1) 



(e[v ]u) ]T 



riGNo 



jn— 1 



nGNo 



L 2 (M 2d ) 



In the following we only estimate the m-th term appearing in the above series: By 
Plancherel's theorem (used for Fe^ri) we obtain 



(Q\y ]u) 



, Sm 



f 2 
Vn 



(3.2) 



L 2 (R 2 



E 



L 2 (R 2ti ) 



nGNo i=0 V ' V ■" 



To proceed further we truncate the above summation in n at an arbitrary N G N. 
We then set I := 2n — j, as well as 



i + (-iy +e 



1 , if j + £ even 
0, else. 



We rewrite (|3.2[) as 



(e»)E 



2iV 



L 2 (R 2d ) 



Ei ! E - ] 2 — " 



3=0 



1 2AT ,j/ 2 

^EthittK^II 



j=0 L 2 



eV^ 



L 2 (R 2 



In the second step we have used Holder's inequality and the following elementary 
result 



J'- 



< 



1 



with [z] € No denoting the biggest natural number smaller or equal to a;. By p.4j) . 
we have that 

sup |^ m (<5Vb)(x,7?)| = 2^ sup |fl^Vo(i/)|. 



Hence, 



(e[v ]u)J2 



n=0 



2^ <$l* 



L 2 (R 2 



< 



E 



d^ 2 



2^ 07 ri+ii i - 2^ 



g the 

1 / 



\k m \ 3 V 



\km\ J V() llii(Kd) 



2£ ,f /2 



Now we note that (by decomposing the remaining sum into its even and odd indices) 
we obtain 



=0 L 2 J ■ j=0 J 



mi I e dfe ™/ 4 . 



We finally arrive at 



(e[v ]u)J2 



n=0 



< 



2e d ' A 



L 2 (R 2 



' i + Vd\k 



™l 1 e dfc ™/ 4 K 



|| e«-w | 



L 2 (R 2d )' 
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for any m = 1, . . . , d and for all N > 1. Thus, (|3.1[) implies 

I ™.) E , ^ » »,,,, E II (l + ^) <** IL, 

nGNo m=l \ / 

which concludes the proof. □ 

Now we are ready to prove our first main result. 

Proof of Theorem [IJ We shall first show assertion (i) : Any solution G TL of 
(II. 6|) that is subject to the normalization JJ dxd£ = 1, satisfies the unique 
decomposition = fi + w* with w* G TL . Therefore, we consider the following 
fixed point iteration for w* : 

T:7Y J -^7Y- L , w n —i i > T(w n —\) = w n , 
where w n G TL solves 

Lw n = XQ[Vo](w n -i + /t). 

To be able to apply Banach's fixed point theorem, we have to prove that the map- 
ping T is a contraction on TL . To this end we write, for any w n —\,w n —\ G TL , 

II w„ - u5 n || w x = || XL^ 1 Q[Vo](w n -i - w n -i) \\h± 

and estimate 

\\w n - w„\\ n ± < |A| || \\ B{n ±j \\Q[V }(w n -i - w n -i) \\ H ± ■ 
From the Propositions 12.41 and 13.21 we consequently obtain 
n ~ ii ^ Co|A| „ ~ 

for all potentials Vq, which satisfy (jl.llj) . Since |A| < a/Co, there exists a unique 
fixed point w* = T(w*) G Thus, the unique (stationary) solution of (jl.6| is 

obtained as Woo = /i + w* , and it holds // R2 d Woo dxd£ = 1. 

Note that the obtained solution Woo G TL C 6>'(R 2d ) satisfies (| 1 .6[) . at least in 
the distributional sense. Moreover, 9[Vb]woo € H and Lwoo G H~ 2 (R 2d ) and thus 
(|1.6p also holds in H~ 2 (M. 2d ). To explore - a-posteriori - the regularity of Woo, we 
rewrite (|1.6[) in the following weak form 

J Jm? d 

for any <p G ff 1 (R 2d ), where 

F{Woo) ■= Woo - divx^Woo) + div^^Woo + 2£«; 0O ) - A6[Vo]ww 

Clearly F^w^) G i?" 1 ^ 2 ^ and thus G ff 1 ^) follows. Moreover, since 
*>cc) £ iL( R2d ). w e also have Wco G fff oc (M 2<i ). 

To show assertion (ii) of Theorem Q] we first note that <9[Vq] is a bounded per- 
turbation of the generator L on TL and thus admits a unique mild solution 
w G C([0,oo),H), c.f. Corollary 12.31 Since 0[Vb] maps H into H -1 , we also know 
that along this solution the mass is conserved, i.e. J J w(t) dxd£ = 1, for all t > 0. 

Next, consider the new unknown v(t) := w(t) — w^. Due to mass conservation 
v(t) G TL , for all i > 0, and we also have 



d t v = Lv - Q[V ]v, 
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since is a stationary solution of (jl.ip . Using this, one easily computes 

4lK*)ll« = 2 // ^-(Lv-Q[V ]v)dxdC 
at JjR2d p 

= 2 {I ^(L s v-e[V ]v)dxd^ 
J Jm.2d p 

where we decompose L according to (|2.2j) . having in mind that L as is anti-symmetric 
in H. Recall now that L s < —a on Ti.- 1 , c.f. Proposition 12.41 Together with 
Proposition 13. 2[ this consequently yields 

^ < 2(C |A| - o-) ||«(t)||^, Vi>0. 

The exponential convergence then follows by a Gronwall estimate and the theorem 
is proved. □ 

Remark 3.3. Due to the mass normalization J J dxd£ — JJ pdxd£ = 1, the 
fixed point w* must take both positive and negative values. Thus, Woo — p + w* 
may, in general, also take negative values. 

Proof of Theorem^ We start with assertion (ii), which follows from the fact that 

Vp G T 2 : ||p - || p(; •) \\ L2 = (2ir) d / 2 \\ w(; •) \\ L2 < K\\ w(; -)\\ n , 
for some K — K(d) > 0. Thus, we infer 

p(t) ^ poo in T 2 

with the exponential rate obtained from Theorem 1 (ii). 

To prove assertion (i) we consider the transient equation (II. lj) as an auxiliary 
problem: Choose any po G if such that tr po = 1 and the corresponding wo G TL. 
Due to the results on the linear Cauchy problem given in [8] we know that 
gives rise to a unique mild solution p G C([0, oo); if ), satisfying trp(i) = 1, for 
all t > 0. The Banach-Alaoglu Theorem then asserts the existence of a sequence 
{AijneN C M+ with t n — > oo, such that 

p(t n ) "-12? p m ^ x weak-* 

for some limiting p G 1\. The already obtained T 2 - convergence of p(t) towards 
Poo € X 2 implies p^ = p G Ti. And the uniqueness of the steady state yields the 
convergences of the whole i-dependent function p(i) — > p in Ti weak-*. Finally, 
we also conclude positivity of the operator p,^ by the T 2 - convergence and the fact 
that we already know from [5]: p(t) > 0, for all t > 0. 

It remains to prove trpoo = 1. To this end, we recall that for any p G the 
corresponding kernel 

&(x,r)) : =p(x + | , 
satisfies G C(K^, L|(M^)), see pQ, and it also holds 

(3.3) trp= / i?(a;,Q)dx. 

Further, note that $(x, 77) = (J 7 ^ n w)(x, rj) = w(x, if), by (|1.2p . On the other hand, 
for any w G 7i we know that u> G C(R^, L 1 (R^)), due to the exponential L 2 -weight 
p -1 in x G M. d and a simple Sobolev imbedding w.r.t. the variable rj G M d . Hence 
the normalization condition dxd£ = 1 implies trp^ = 1, via (|3.3|) . and 

assertion (ii) is proved. 

Finally, we prove claim (iii) by first noting that the T 2 -convergence of p(t) implies 
convergence in the strong operator topology. Thus, having in mind that ||p(f)|| T = 
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|||0o || = 1) we infer from Griimm's theorem [TH] that p(t) also converges in the 
Ti-norm towards poo. In summary, this concludes the proof of Theorem [H □ 
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